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SUMMARY 


A single-wire magnetic leader cable laid along the path the aircraft is to 
follow is a potential source of guidance information during rollout, turnoff, and 
taxi. This report contains a theoretical analysis of the errors produced in such a 
system by the geometry of the cable layout and by a nonzero aircraft attitude. It 
was found that neglecting the aircraft attitude in processing the leader cable sig- 
nals to obtain the aircraft heading can cause significant errors, even for small 
attitude angles. Displacement from the cable, on the other hand, can be obtained 
with reasonable accuracy without compensating for aircraft attitude. 

Analysis of a rectangular loop geometry for the cable shows that for practical 
loop widths, a nonlinear calibration curve is required to obtain accurate estimates 
of aircraft heading and lateral displacement relative to the cable, although a linear 
approximation may be sufficiently accurate for many applications. 

From analysis of a practical runway turnoff geometry, it was found that the 
errors produced by the turnoff are small for the estimates of displacement and are 
largest for the estimates of heading at the beginning and end of each turn. The 
magnitude of the error was not particularly sensitive to the turnoff angle, but the 
heading error did increase with a decreasing turn radius. In all cases, the errors 
were minimum at or near the cable. 

Although the leader cable concept was originally formulated for the hypothetical 
case of an infinite straight wire, results of this analysis show that the concept 
remains viable for a practical closed-loop runway-turnoff-taxiway configuration. Of 
course, the effect of all leader cable sensor errors on total guidance and control 
system performance must be evaluated via computer simulation and flight tests. 


INTRODUCTION 

The new microwave landing system (MLS) being developed by the Federal Aviation 
Administration is designed to provide sufficiently accurate guidance during final 
approach and touchdown to permit aircraft to land under conditions of severely 
reduced visibility. In addition, the azimuth and DME functions of the MLS will pro- 
vide some guidance during rollout under these same poor visibility conditions. How- 
ever, guidance during turnoff and taxi is currently available only through visual 
means and thus is inadequate during periods of extremely low visibility. Further- 
more, the MLS, together with advanced avionics and air-traffic-control techniques, 
has the potential to increase the runway capacity by decreasing the spacing between 
landing aircraft. Decreased landing intervals demand that runway occupancy time be 
decreased, that is, that landing aircraft be cleared from the runway rapidly. One 
technique for reducing runway occupancy time is the use of high-speed turnoffs. 

Safe and efficient use of high-speed turnoffs requires the availability of turn- 
off guidance, particularly for use with automatic controls. Clearly, if the MLS and 
high-speed turnoffs are to be used to their peak potential, appropriate guidance 
information during turnoff and taxi must be made available. One potential source of 
such guidance is the magnetic leader cable. The cable, or wire, would be buried in 
the runway, turnoff, and taxiway along the path that the aircraft is to follow. (See 



fig. 1 . ) A sensor in the aircraft would detect signals from the cable and provide 
guidance to the pilot or automatic control system. 



Figure 1.- Rollout and turnoff guidance using a magnetic leader cable. 


As early as the 1950 ’S/ the magnetic leader cable was investigated by the 
British (refs. 1 and 2) as a source of guidance during final approach and landing, 
but it was never adopted as an operational system. One of the factors limiting its 
usefulness in that application was the difficulty in obtaining accurate guidance 
during final approach prior to touchdown. This difficulty would be avoided in the 
present concept by utilizing the system only after touchdown. Other applications 
considered for the magnetic leader cable include guidance of automobiles (ref. 3) and 
special-purpose military vehicles (ref. 4). 

Work on the magnetic leader cable for aircraft guidance has been revived as part 
of the NASA Terminal Configured Vehicle (TCV) Program. One facet of this effort 
involves the investigation via computer simulation of aircraft performance during 
rollout and turnoff using MLS and magnetic leader-cable sensors. Early portions of 
this work have been reported in references 5, 6, and 7. A second facet is the devel- 
opment of a magnetic leader cable sensor system. Although two-wire cable configura- 
tions have been considered for other applications (ref. 3), a single-wire system was 
selected for the current application to simplify the task of embedding the cable in 
the runway and taxiway along the desired path. One part of this development is a 
theoretical analysis of a practically configured cable and sensor system. Such a 
theoretical analysis is the subject of this report. Previous analyses have assumed 
an infinitely long, straight wire; the present analysis investigates the effects of 
aircraft attitude, a finite loop size, and a typical turnoff geometry on the poten- 
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tial accuracy of the system. Together with system performance requirements and other 
data, the results of this analysis can then be used to determine such things as mini- 
mum loop width and the need for aircraft attitude compensation. 


SYMBOLS 


2 

A area of coil, m 

a-j,a 2 /a^,a^ intermediate variables used in computation of magnetic intensities 

^4' respectively (see eqs. (21), (23), (25), 

and (27) ) , m^ 

->• 2 
B(t) magnetic induction field, Wb/m 

Cartesian components of B, Wb/m"^ 

d perpendicular displacement from cable in xy-plane, m 

a incremental length of current element, m 

d^ ,d 2 estimate of d defined by equations (18a) and (18b), same as y^ and y 2 
before turnoff 


LL (ijNAI 

H magnetic intensity. A/m 

Cartesian components of H, A/m 

X y z 

magnetic intensity produced by current in sides 1, 2, 3, and 4, 
respectively, of rectangular loop. A/m 

I peak current in cable, A 

I(t) time-varying current in cable, A 


i ,i ,i unit vectors along Cartesian axes 

X y z 

A A ^ 

i^,i 2 ,i 3 unit vectors perpendicular to planes of coils 1, 2, and 3, respectively 

-► 

i unit vectors in 0-direction in cylindrical coordinates 

0 

L length of cable loop in x-direction, m 

Lg length of straight segment of cable between two turns of turnoff, m 

L^ x-coordinate at start of turnoff, m 

N number of turns in coil 

P point at which magnetic field is calculated 

R^,R 2 radii of first and second turns, respectively, in turnoff, m 
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« >1 


r radius vector from origin to point P, m 

vector from current element I dl to point P, m 
S surface area, 

distance between runway and taxiway, m 

s curvilinear distance along path parallel to and at a distance d from 

cable , m 

t time, sec 

peak voltage output of coils 1, 2, and 3, respectively, V 
peak voltage at sensor's y^ output, V 

voltage equivalent to Vy^ in x' ,y* ,z* coordinate system 


v^(t) 

W 

x,y,z 
x' ,y’ ,z' 




yi'^2 

A 

^1 


0 , (|), <\) 


time-varying voltage output of coil 1 , V 

width of cable loop in y-direction, m 

Cartesian coordinates of point P or of coils 

coordinates where x* = -x, y' = -y, and z* = z 
Cartesian coordinates of cable current element 

estimates of coil displacement from cable defined by equations (17a) 
and ( 17b) , m 

estimates of coil displacement from cable defined by equations (18a) 
and ( 18b) , m 

estimate of coil displacement from cable using nonlinear calibration 
curve , m 

A 

A 

estimate equivalent to y-| in x*,y',z* coordinate system 

angles of pitch, roll, and yaw, respectively (9 also used as angular 
coordinate in cylindrical coordinates), rad 


(i)^(t) 


turnoff, or exit, angle, rad or deg 

-7 

permeability of free space, 4it x 10 Wb/A-m 
radial coordinate in cylindrical coordinates, m 
magnetic flux linking coil 1, Wb 
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(\> aircraft heading relative to cable (equivalent to yaw during rollout ) , 

rad or deg 

^ estimate of c|; defined by actuation (17d)/ rad or deg 

A 

4> estimate of <|; defined by equation (18c)/ rad or deg 

0 ) angular frequency of current in cable/ rad/sec 

Abbreviations : 

DME distance measuring equipment 

MLS microwave landing system 

TCV terminal configured vehicle 


DEVELOPMENT OF THE THEORY 

Consider an infinitely long/ straight conductor coincident with the x-axis 
carrying a current I(t)/ as depicted in figure 2. The magnetic leader cable uti- 
lizes a sinusoidal current in the audio frequency range/ say 100 Hz to 2500 Hz* 
Because of the extremely long wavelengths (greater than 100 km) and short distances 
(less than 5 km) involved/ the propagation time can be neglected. Thus in order to 
simplify the mathematics/ the time dependence of the current and of the resulting 
magnetic field will be suppressed/ and the quasi-static approximation will be used in 
the magnetic field calculations in the remainder of this report. 


z 



/ 


Figure 2.- Geometry for infinite straight cable. 

Returning now to figure 2/ the magnetic intensity H at point P can be com- 
puted using Ampere's circuital law as follows: 
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( 1 ) 


f H • dt = I 
C 

Because the field is symmetrical about the conductor and is everywhere in the 
©^direction, equation (1) simplifies to 


-► 


H = 


I 

2np 


( 2 ) 


In Cartesian coordinates , equation (2) becomes 


= 0 


__ I z 

”y 2it 2 ^ 2 

^ y + z 


H 

z 


I 

2n 




2 

z 


(3a) 

(3b) 

(3c) 


Suppose there were three ^mutually orthogonal coils located at point P such 
that the axes i^z ± 2 ' ^3 coils were parallel to the coordinate axis 

Xz yz and z, respectively- From Faraday's law the voltage induced in coil 1 by the 
magnetic field is 


d(|) (t) 

^ 

If coil 1 has number of turns N with area A and the cable current is sinusoidal 
with angular frequency tOz the flux linkage is given by the integral 


4)^(t) = N / ]l(t) ♦ i-j dS = NAB^ cos o)t 


(5) 


The voltage induced in coil 1 is then 


v^(t) = oiTAB^ sin ojt 


or in terms of peak values 


= 0>NAB = OdNAH 
1 X ^O X 


Combining equations (3a) and (6b) z we see that 


(6a) 


(6b) 


(7a) 
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Siitiilarly for coils 2 and 3, we find that 


ix cjnai 

V ^ 

2 2n 2^2 

y + z 

and 


(7b) 


V 


3 


U OJNAI 
^o y 

2n 2^2 
y + z 


(7c) 


Now take the ratio of voltages to V 2 and solve for 


y. Thus, 


y = 



( 8 ) 


If the coil height z is known, then the displacement y of the coils from the 
magnetic leader cable can be determined by measuring the coil voltages V 2 and 
and using equation ( 8 )« Note that in this case the displacement y is a linear 
function of the ratio of voltages ^ 3 /^ 2 * Th® displacement y can also be deter- 
mined via an alternative technique which rec[uires a knowledge of the current I in 
the cable but which does not require that the height z be known* This can be seen 
by finding the ratio of to V 2 ^ + as follows: 


V. 


2 2 
2 3 


27ry 

LL toINA 
^o 


(9) 


Equation (9) can be solved for 


y to obtain the following: 


y 


\i COINA V. 

o 3 


2tz 



+ V. 


( 10 ) 


Note that this e 3 q>ression for y 


is a linear function of the ratio 



V3^) 


Suppose now that the coils are rotated in the xy-plane such that i^ makes an 
angle (\> with the x-axis, as in figure 3* Now the voltage output of the z-coil 

remains the same, but (since equals zero) the voltages and V 2 become 


V 


1 


cdNAB^ 


sin <\> 


( 11 a) 
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V- = aWAB cos (p 

2 y 


(11b) 


z 



Taking the ratio V^/V 2 , we find that 


V 


1 


V 


2 


tan (p 


Equation (12) can be solved for cp to give 


( 12 ) 


cp = arc tan (V^/V 2 ^ 


(13a) 


which, for small values of cp, can be approximated by 


cp 



(13b) 


Note that in this approximation cp is a linear function of the ratio V^/V 2 * 
Equation (8) for the displacement now becomes 


zV^ cos cp 



However, if cp is small, equation (14) can be approximated by equation (8). 

We have shown conceptually that for an infinitely long, straight cable carrying 
a current 1, the lateral displacement y (relative to the cable) of a set of three 
mutually orthogonal coils can be determined by measuring the voltage output of the 
appropriate coils and processing these measurements with either of two algorithms* 
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The use of one algorithm requires that the height of the coils be known, while use of 
the other algorithm requires that the magnitude of the current in the cable be known. 
Furthermore, the angle of rotation of the coil in the xy-plane can be determined by 
using the coil voltages and processing them with a third algorithm. Conceptually it 
follows that if the coils are rigidly attached to an aircraft and the magnetic cable 
is laid on the runway center line, the sensor can be used to measure the displacement 
of the aircraft from the center line and the heading of the aircraft relative to the 
center line. In the remainder of this report we will examine the effects on system 
accuracy when the linear relationships developed above are used to process signals in 
practical situations involving nonzero aircraft attitudes, rectangular loops, and 
different turnoff geometries. 


EFFECTS OF AIRCRAFT ATTITUDE 


The previous development has assumed that the z-coil has its axis in the 
z-direction. With the coils in an aircraft this may not always be true, even when 
the aircraft is taxiing. Suppose the coils have been rotated through a yaw angle 
([;, a pitch angle 0, and a roll angle (|), in that order. The voltage outputs of the 
coils now become 



cos 9 cos (\) 


oNA 


-cos 4) sin ([> 

+ sin <() sin 0 cos <\) 


sin 4> sin <\) 

+ cos (j) sin 0 cos (jj 


cos 0 sin 4; 


cos (J) cos (\f 

t sin <|) sin 0 sin 


-sin ^ cos (\) 

+ cos (t) sin 0 sin (\) 


-sin 0 


— — 
B 

X 

sin (() cos 0 


B 

y 

cos (}) cos 0 


B 

z 



^ - 


(15a) 


For small angles, equation (15a) can be approximated by 




r n 


— — 



e- 

1 

CD 


B 

1 




X 

< 

to 

^ ONA 

-4, 1 <i) 


B 

y 

^3 


CD 


B 

z 

^ — 


-J 


— — 


Combining equation (15b) with equations (3), the coil voltages are 


(15b) 


— — 


r n 


— 



CD 

I 

•a- 


0 





-z 

^2 

= G 

-cP 1 4) 


2 ^ 2 





y + z 

^3 


e - 4 . 1 


y 

2^2 

_ 


■— 


y + z 


(16) 


9 



where G = 


\x oNAI 
2n • 


Equation (16) can be solved for the desired variables y and <\)i 



G(i + 

Y2 = 2 2 " 

v/ + V32 


(|)Z 




GV, 


+ V, 


- (})Z 


cp = 



(17a) 


(17b) 


(17c) 


The solution for ^ requires a knowledge of y, which is an unknown to be deter- 
mined* In order to compute use for y the value y^ calculated in equation 
(17a) - or alternatively y 2 • Then 



These algorithms (eqs. (17a), (17b), and (17d)), which contain small angle approxima- 
tions, can be used to compute estimates of the displacement and heading if the air- 
craft attitude angles 0 and ^ are available. 

During rollout, turnoff, and taxi, the aircraft attitude angles should be close 
to zero. Extreme values for yaw, pitch, and roll probably will not exceed ±10®, ±5°, 
and ±5®, respectively. Within these limits for the attitude angles, what then are 
the errors in the estimated displacement and heading introduced by the small-angle 
approximations in equations (15b) through (17d)? These errors have been computed for 
several combinations of attitude angles for a coil height of 2.5 m; the errors in 
y^ and ^ are plotted in figure 4. From figure 4(a) it can be seen that for the 
attitudes used in computing these results, the errors in y^ produced by the small- 
angle approximations never exceed 0.39 m for displacements y up to ±30 m. The 
error in y 2 never exceeds 0.05 m and is not plotted. The effects of these errors 
on the performance of an automatic rollout and turnoff system must be evaluated via 
simulation and test. However, it is expected that these errors can easily be toler- 
ated. The error in ^ produced by the small- angle approximation is seen from fig- 
ure 4(b) to be less than 0.6® for all the cases considered and is less than 0.01® 
for (\) = 9 = (() = 1®. Thus, it would seem that the error in the estimate ^ may be 
tolerable. 
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(a) Error in estimate . 
Figure 4.- Error in estimates 




0, 

<P, 


deg 

deg 

deg 

o 

10 

5 

5 

□ 

5 

5 

5 

o 

1 

1 

1 

A 

-5 

-5 

-5 


y^ and (JJ due 



(b) Error in estimate 
to small-angle approximations. 


The complexity in instrumentation and computation would be reduced considerably 
if the aircraft attitude could be ignored/ that iS/ if the displacement and heading 
could be calculated from the basic algorithms (eqs. (8)/ (10)/ and (13b)) as follows: 



^2 


^LqCONAI 

2n 




2 2 
^2 + ^ 3 ^ 


4> = 


Xi 


(18a) 


(18b) 


(18c) 


The errors in the estimates computed with these algorithms are plotted for several 
values of attitude in figure 5. From figure 5(a) the error in y-^ is much larger 
than the error in y^ / but the error does decrease dramatically nearer the cable 
(as |y| decreases). From figure 5(b) the error in Y 2 is much smaller than the 
error in y-j . (Note the difference in scales in the two figures.) It would appear 
that the y 2 data are potentially useful in an automatic system. From figure 5(c) 
the error in c|; is very large. In fact/ $ is probably unusable because the errors 
are large enough in some cases to produce an incorrect algebraic sign for a con- 

dition which would produce an unstable guidance system. 
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(a) Error in estimate . 
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(b) Error in estimate y 2 « 
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y\ ^ 

Figure 5.- Error in estimates y-j , y 2 ^ and 4; due to neglecting aircraft attitude. 


EFFECTS OF CLOSED-LOOP GEOMETRY 

The previous equations were developed assuming an infinitely long, straight 
cable. Obviously, a practical installation will not meet these assumptions- The 
cable will not be infinite and may not be straight, and a return path for the current 
must be provided. 

Consider a cable laid along the runway center line with the return in the form 
of a closed rectangular loop, as depicted in figure 6. The magnetic intensity pro- 
duced at point P by a current I in this lo<^ can be found by using the Biot- 
Savart law to find the incremental intensity dH due to a current element I dl 
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and then performing a line integration around the closed loop to get the total inten- 
sity H; that is. 



Figure 6.- Geometry of rectangular-loop cable. 


The integral in equation (19) can be evaluated in closed form by expressing the 
integrand in rectangular coordinates and integrating separately over each of the four 
sides of the rectangle. For side 1 the results are 


^L/2 

*^V2 [(X - + y2 + 


^2j3/2 


4na 


+ yiz) 


where 


X + 


X - ^ 

^ 2 




X <-|)2 + a ,: 


i 


X - 1)^ . a,^ 


( 20 ) 


2 2 2 
= y + z 


( 21 ) 
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similarly, for sides 2 


3, and 4: 


^2 

4ita, 


zi - (y - W)i 


A 


X + § 


^(x + i)2 + a2^ 


^ _2_ 


ilF« 


2 2 


( 22 ) 


where 


2 2 2 
= (y - m ^ 


(23) 


S, - * (x f i)ll 


y - W 




2 2 
W)"-H a3^ 


(24) 


where 


.(x.f) 


2^2 
+ z 


(25) 


-*■ 

«4 = 


4ii;a^ 


- (* - 


y - W 


^(y - 


2 2 
W) + ^4 


(26) 


where 


- (» - 1 )= 


\2 2 
= I X - r ) + z 


(27) 


The total magnetic intensity is given by the sum of the four components: 


«2 «3 ^ «4 


(28) 


Suppose that three orthogonal coils are used to measure the components of the mag- 
netic field produced by the rectangular loop. Using the basic algorithms in equa- 
tions (18),^let the output voltages V-j , V 2 , and be used to compute the 

estimates y-| , y 2 r and These estimates will be in error because equations (18) 

were developed assuming that the magnetic field was produced by current in an infi- 
nite, straight wire rather than by a closed loop. Let us next examine the magnitudes 
of these errors as functions of the various parameters. 


Variation in Error With y-Position 

Figure 7 contains plots of the errors in y^ and y 2 caused by a rectangular 
loop 3000 m long and 100 m wid6. The errors are plotted as a function of y for 
two values of coil height, z = 1 and 6 m. These heights should be near the extremes 
of heights for coils mounted in the nose of current aircraft. These data are for 
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Displacement from cable, y, m 


A A 

Figure 7.- Error in estimates and Y 2 for a 100-m-wide rectangular loop. 


A 

X = 0; that is, the coil is equidistant from the ends of the loop. The error in y-j 
is larger for positive than for negative values of y because the return side of the 
cable loop is on the positive side (at y = 100 m) . There is very little difference 
in this case in the error for the two different heights. The error in y 2 is more 
nearly symmetrical about y = 0 than the error in y^ and is of the opposite sign. 

It is worth noting that the error is smallest near the cable (y = 0) . This is an 

important factor since the pilot or the automatic control system will be trying to 
follow the cable. The error in ^ in this case is zero because the coils are 
located at x = 0. For all the data in figure 7, the true relative heading is zero 

((|; = 0®). Results computed for (f) - 5^ show errors very close to those in figure 7. 

The errors in figure 7 are bias-type errors which cannot be reduced by filter- 

ing. However, these errors could be reduced, or even eliminated, by using a nonlin- 
ear calibration curve to convert the sensor output voltages to position and heading 
estimates instead of using the linear relationships from equations (18). Of course, 
such nonlinear relationships increase the complexity of the system. The calibration 

curve which accurately relates the sensor output V = zV /V to the estimate yi 

y ^ 3 2 

is shown in figure 8 for a loop width of 100 m. This calibration will be less accu- 
rate as the loop width deviates from 100 m. 
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Figure 8«- Nonlinear calibration curves for estimating displacement y for 

a 100-m-wide rectangular loop. 


There is an additional problem to consider. The calibration labeled 


figure 8 and the error curves in figure 7 are for the return wire located at 
y = 100 m and for an aircraft landing in the negative x-direction. Consider an 
aircraft on the same runway landing in the positive x-direction. For this case, 
define a new set of coordinates x',y*,z*, where x* = ~x, y* = -y, and z’ = z. 

Now the return wire is located at y’ = -100 m, and it can be shown that the calibra- 
tion curve relating the sensor output V' to y* is the curve V,, reflected 

1 y-i 


about the origin, as shown in figure 8. Thus, to accurately remove the loop effects 
with a nonlinear calibration curve requires that the sensor store two curves and that 
the proper curve be selected for the given landing situation. Alternatively, the 
sensor could compute the reflection about the origin of a single calibration curve. 
However, it should be remembered that such accuracy may not be required by an auto- 
matic rollout control system. The accuracy requirements must yet be determined via 
computer simulation and flight testing. 
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Variation in Error With Loop Width 


The errors could also be reduced by increasing the width of the loop and 
thereby decreasing the effect of the current in the return path. Figure 9 shows the 



Displacement from cable, y, m 


Figure 9.- Error in estimates y^ and ^2 ^ 300-m-wide rectangular loop. 


same type of errors as figure 7 except that the loop width is 300 m instead of 
100 m. As expected/ the errors for the 300-m loop are smaller. For example/ at 
y = 20 m the error in the estimate y-j is 7 m for the 100-m loop and less than 2 m 
for the 300-m loop. The errors in y-j and Y2 plotted in figure 10 as a func- 

tion of loop width W; it can be seen that these errors decrease very rapidly with 
increasing width for W < 200 m. For large values of W (above 600), the rate of 
decrease in the error becomes very small. The obvious problem in increasing the loop 
width at a typical airport is the problem of having a place to lay the loop. 


Variation in Error With x-Position 

The results in figures 7 through 10 are for x = 0, and the errors are due 
primarily to the effects of the return side of the rectangular loop. Near the ends 
of the loop, errors are also produced by the current in the end. In figure 11 the 
errors are plotted as a function of the distance along the cable x for two values 
of Yf y = 5 and 20 m. The loop is 3000 m by 100 m, and the coil is at a height 
of 2.5 m. From a value of 7 m at x = 0, the error in y-j at y = 20 m 
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error 



Loop width, W, m 

Figure 10.- Error in estimates and ^2 function of loop width. 

y = 20 m. 


increases to 11 m at a distance 50 m from either end. The effect of the ends on (|; 
is even more pronounced. From a value close to zero (less than 0.1®) along most of 
the cable, the magnitude of the error in c|^ increased rapidly within a few hundred 
meters of the end to a value of 2.5® at 50 m from the end. Unlike the error as a 
function of y, these end-effect errors cannot be readily reduced by using a nonlin- 
ear calibration curve. To reduce end effects, the cable must be laid such that the 
ends are some 200 to 300 m beyond the active portion of the cable, that is, that part 
of the cable which the aircraft will follow. 






EFFECTS OF TURNOFF GEOMETRY 


The rectangular loop is configured to provide guidance in a straight line along 
one side of the rectangle, such as along the runway center line during rollout. How- 
ever, it is also desired to provide guidance during turnoff and taxiing. Let us next 
examine the errors introduced by laying the cable in a runway, turnoff, and taxi way 
conf i guration . 

Consider the geometry in figure 12. The magnetic intensity at point P is 
found via line integration of the Biot-Savart law, as in the case of the rectangular 
loop. The contribution to the total field of the current segments 1 through 5 has 
the same form as for the rectangular loop. The integrals describing the fields due 
to the linear and circular arc portions of the turnoff (segments 6, 7, and 8) were 
not solved in closed form, but rather were evaluated by numerical integration. 



X 


Figure 12.- Loop geometry for turnoff. 



These fields were then used to compute^ the voltages in the coils, which in turn were 
used to compute the heading estimate and the displacement estimates d^ and 

<± 2 , using the basic algorithms in equations (18). The variables d^ and d 2 are 
estimates of d, the linear displacement from the cable in a direction perpendicular 
to the cable. The variables d, d^, and d 2 correspond to the previously used y, 
y^ , and Y 2 # but the change in nomenclature is necessitated by the fact that the 
displacement is not always in the y-direction during turnoff. 

The dimensions of a typical runway- turnoff -taxiway configuration are listed in 
the following table: 


Loop length, L, 3000 

Loop width, W, m 100 

Runway- taxi way separation, ^ep' * 150 

Start of turnoff, L^, m 150 

Turnoff angle, 6^, deg 30 

Radius of first turn, R-j , 550 

Radius of second turn, R 2 , m 250 

Distance between turns, L^, m 85.6 

Height of coils, z, m 2.5 


Unless otherwise noted, these dimensions were used in the computations that follow. 


Variation in Error With Along-Track Position 

The errors produced by the cable geometry were computed as the coils were moved 
along a path parallel to and at a distance d away from the cable. These errors are 
plotted as functions of s in figure 13 for values of d = 5 and 20 m. The 



(a) Error in estimates d^ and d 2 . 
Figure 13.- Error due to turnoff geometry. 
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(b) Error in estimate ([>. 


Figure 13.- Concluded. 


variable s is the curvilinear distance along the path starting at x = L/2 and 
moving in a negative x-direction; it is analogous to the variable x in the error 
plots for the rectangular loop. 

The following observations can be made about these errors: 

1. The errors decrease significantly as the aircraft (coils) approaches the 
center line (the cable ) , just as in the case of the rectangular loop cable geometry. 

A yv 

2. The errors in the displacement estimates d^ and ^2 sltq close to their 
values for a rectangular loop in areas away from the turns. In the region 

1855 < s < 3000 the cable is 100 m from the return wire, and the errors in figure 13 
are nearly the same as the errors in figure 11 for a rectangular loop 100 m wide. In 
the region 0 < s < 1350 the cable is 250 ra from the return wire, and the errors in 
this region are nearly the same as the errors for a rectangular loop 250 m wide. The 
error in one region, but not both regions, can be removed by a single calibration 
curve. Conceptually, multiple calibration curves could be used to reduce the errors 
in both regions at the expense of increased complexity. Of course, the difference in 
the errors between the two regions could be reduced by increasing the width of the 
loop. 


3. The errors in the displacement estimates d-j and 6.2 smoothly in the 

areas of the turns. 

4. In the plot of the error in the estimate <|; (fig. 13(b)), spikes occur at 
the beginning and end of each turn. The magnitudes of these spikes are larger for 
the larger values of d. 

A 

5. In the regions away from the turns, the error in <|; is nearly zero, as in 
the case of the rectangular loop. It is worth noting that the sign of the error in 
4> at the beginning and end of the turns is such as to provide some lead, or 
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predictive r information going into and coming out of the turns. The usefulness of 
this information has not been evaluated. 


Variation of Error With Cross-Track Position 

A /V ^ 

In figures 14 and 15 the errors in the estimates d^ , 6 ^ 2 ^ and (\> are plotted 

as functions of d. In figure 14 the path is located at the beginning of the first 
turn, and in figure 15 the path is located at the beginning of the second turn. 



Figure 14.- Error in estimates d^ , 


d^, and <\> at start of first turn. 


Figure 15.- 
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Error in estimates d^ , 


d 2 f and 


(j; at start of second turn. 



I 


These locations correspond approximately to the locations of the spikes in the 
4r-error in figure 13(b). We observe that the magnitudes of the errors generally 
decrease smoothly with decreasing magnitude of d and are minimum near d = 0. One 
exception is the local maximum for the magnitude of the (|/-error near d = 20 m in 
figure 15. 

Variation of Error With Turnoff Angle 

Consider now the effect of the turnoff, or exit, angle 9^ on the magnitude of 
the errors during turnoff. In figure 16 the maximum errors in^the estimates during 



(a) Error in estimates d^ and d^- 

Figure 16.- Maximum error magnitude during turnoff caused by cable geometry 

as function of turnoff angle. 
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Turnoff angle, deg 


(b) Error in estimate (\)* 
Figure 16.- Concluded. 


the first turn are plotted as functions of the exit angle for d = 0, 5, and 20 m. 
The radius of the first turn is 550 m, as used previously. From these graphs 
(figs. 16(a) and 16(b)) it can be seen that for positions close to the cable (small 
values of d)/ the peak errors are not significantly affected by the exit angle for 


R-| = 550 m. At d = 20 m the magnitude of the errors in d-^ and 6.2 increases 
slightly with 0g, but the effect is not large enough to be a determining factor in 
selecting the exit angle. Larger exit angles are not considered in figure 16 because 
angles greater than 35.66® are not possible with the given geometry (R^ = 550 m. 


250 m, and 


= 150 m) 
eo 





Variation of Error With Turnoff Radius 


We now consider the effect of the turnoff radius on the errors in the esti- 
mates. In figure 17 the maximum error magnitudes are plotted as functions of the 
turnoff radius for an exit angle of 30®. Maximum error magnitude is meant to be 
the largest error magnitude encountered as a function of s in the vicinity of the 
turnoff. From figure 17(a), it can be seen that the maximum errors in the estimates 



200 300 400 500 600 700 800 900 

Turnoff radius, , m 


(a) Error in estimates d-^ and d 2 . 

Figure 17.- Maximum error magnitude during turnoff caused by cable geometry 

as function of turnoff radius. 
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Turnoff radius, m 


(b) Error in estimate 
Figure 17.- Concluded. 


d-| and 6.2 are weak functions of the radius R-| for small values of d. On the 
other hand/ the peak error in the heading estimate ^ decreases with increasing 
radius, particularly for larger values of d. To minimize the error in ([>, there- 
fore, R-| should be made large- Of course, other factors, such as acceleration on 
the aircraft and available real estate, affect the choice of R-j . 

SUMMARY OF RESULTS 

This report has examined via a theoretical analysis the errors in a single- wire 
magnetic leader cable sensor system caused by aircraft attitude and by cable loop 
geometry. A set of basic algorithms was developed which could be used to compute 


y and <|; (where y is the aircraft lateral displacement from the cable and c(; is 
the heading relative to the cable) . These algorithms assumed that the cable was con- 
figured as an infinitely long, straight wire; also the equations neglected aircraft 
pitch and roll. 

For the estimate ^2 of the displacement y, it has been shown that the error 
caused by neglecting aircraft attitude may be small enough to be neglected. However, 
for the estimate cp of aircraft relative heading (p, the error caused by aircraft 
attitude is large enough that the heading estimate may not be usable. The attitude- 
induced errors can be significantly reduced by employing the algorithms which include 
the aircraft attitude angles in the computations using small-angle approximations. 

The errors produced by using the basic algorithms when the cable was configured 
as a closed rectangular loop rather than an infinite straight wire were examined. It 
was found that these errors decreased in magnitude with decreasing y and reached 
minimum at or near the cable. The errors increased sharply as the loop width 
decreased below 200 m; but, if required, these errors could be removed for any given 
width by using a nonlinear calibration curve for the computation - in other words, by 
using a more complex algorithm. This calibration curve must be reflected about the 
origin when heading in the opposite direction on the runway. Also, the errors 
increased near the ends of the loop, necessitating that the ends be located 200 to 
300 m beyond the active, or guidance, portion of the cable. It was found that the 
turnoff produced additional errors in the estimates, and these errors also decreased 
in magnitude nearer the cable. Particularly noticeable are the error spikes in the 
heading estimate at the beginning and ending of the turns. It was found that the 
maximum magnitudes of the errors in the turnoff were not significantly affected by 
the turnoff angle for angles in the range of 10° to 35°. The maximimi error magni- 
tudes in the estimates of displacement were not strongly affected by the radius of 
the turnoff for radii between 300 and 800 m. The error in cp, however, did decrease 
significantly with increasing radius for large values of displacement. Although the 
leader cable concept was originally formulated for the hypothetical case of an 
infinite straight wire, results of this analysis show that the concept remains viable 
for a practical closed- loop runway- turnoff -taxiway configuration. 

It should be remembered that there are other sources of error, such as ground 
effects, aircraft metal effects, and sensor noise, which have not been considered in 
this report. Furthermore, the effect of all errors on aircraft performance during 
rollout and turnoff, including those analyzed in this report, must be examined via 
simulation and experimental testing. The sensor performance requirements can then be 
specified. From these requirements, from the test data, and from the results of the 
current analysis, the need for aircraft attitude compensation and for nonlinear cali- 
bration in the sensor processor can be determined. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 5, 1982 
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